Symmetry algebra of discrete KdV equations and corresponding 
differential-difference equations of Volterra type 



Pavlos Xenitidis 
School of Mathematics, University of Leeds, UK 
P.Xenitidis @leeds. ac.uk 



Abstract 

A sequence of canonical conservation laws for all the Adler-Bobenko-Suris equations is derived and is em- 
ployed in the construction of a hierarchy of master symmetries for equations H1-H3, Q1-Q3. For the discrete 
potential and Schwarzian KdV equations it is shown that their local generalized symmetries and non-local mas- 
ter symmetries in each lattice direction form centerless Virasoro type algebras. In particular, for the discrete 
potential KdV, the structure of its symmetry algebra is explicitly given. Interpreting the hierarchies of symme- 
tries of equations H1-H3, Q1-Q3 as differential-difference equations of Yamilov's discretization of Krichever- 
Novikov equation, corresponding hierarchies of isospectral and non-isospectral zero curvature representations 
are derived for all of them. 

1 Introduction 

Infinite hierarchies of symmetries and conservation laws are among the characteristics which most of the well 
known integrable systems possess |[T6l l22l l30l l36l . Their existence, serving as a definition of integrability, has 
been used in the classification of integrable differential equations, eg. H [16l [30l, and differential-difference 
equations ||4l[T3l|27l[34l[35'|, while this approach was recently extended to the case of difference equations Ill4l . 
An alternative approach to this notion of integrability is the existence of either a recursion operator [22] or master 
symmetries (111221, which can be used recursively for the construction of symmetries and conservation laws. 

The above integrabihty criteria can be applied to difference equations as well, for which there exist also some 
other, alternative criteria to establish integrability flOl. Among them can be found the multidimensional consis- 
tency L6j J9„20J and the algebraic entropy [5 |. Adler, Bobenko and Suris classified all the affine linear quadrilateral 
equations which possess the former property |j2] |3], while it was shown in \29] that the algebraic entropy of any 
equation in this class vanishes. Moreover, the symmetry analysis of these equations 1241 l28l revealed that one 
of their symmetries plays the role of master symmetry with the help of which a hierarchy of symmetries was ex- 
plicitly constructed in [31]. Moreover, a recursion operator was recently found for the most general equation in 
this class along with some canonical conservation laws [[TTlfTSl . while some other conservation laws were found 
constructively in Il23l and by implementation of the Gardner method in ||25| . In this paper, we focus on the Adler- 
Bobenko-Suris (ABS) equations and explore further their integrability features. Among these integrability features 
are an infinite hierarchy of canonical conservation laws and a hierarchy of master symmetries for equations H1-H3 
and Q1-Q3. 

More precisely, first we present a derivation of a hierarchy of canonical conservation laws for all the ABS 
equations. This derivation is similar to the approach used in |T| and it is based on the use of the first canonical 
conservation law, found in [17], and the first master symmetry, given in [31]. In particular, we introduce a potential 
in terms of which the conservation law can be given as the compatibility condition of a system of two, first order 
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difference equations. Extending the action of the master symmetry to this potential results to the introduction of 
a new potential via another system of first order difference equations. The compatibility condition of the latter 
system is the next, higher order, canonical conservation law of the original ABS equation. Repeating the same 
argument with the new system, one is led to the next canonical conservation law, and, successively, a hierarchy of 
canonical conservation laws can be constructed. 

The introduction of these potentials allows us to write in a simple and elegant form the higher order master 
symmetries of the ABS equations H1-H3 and Q1-Q3. More precisely, while the first master symmetry is local, the 
second member of this chain, and subsequently all the higher order members, depend explicitly on the potentials 
of canonical conservation laws. In spite of the fact that all the higher order master symmetries are nonlocal, it 
is proven that their commutators with the local generalized symmetries result to local generalized symmetries. 
This is in accordance with the continuous case [26 ] where sufficient conditions guarantee the locality of these 
commutators. Actually, there exist obvious similarities in the form of the master symmetries constructed here and 
their forms given in Proposition 2 of [26 1. 

For two of the ABS equations, namely equations HI and QIq, we find the algebra of symmetries and master 
symmetries in each lattice direction. Specifically, we find explicit expressions for the characteristics of all sym- 
metries, as well as how they act on the potentials of canonical conservation laws. Subsequently, we compute the 
commutators among the members of each hierarchy and find that they form centerless Virasoro type algebras. 
In particular for HI, we are able to find the commutators among the two sub-algebras and, thus, to describe its 
symmetry algebra. 

From the interpretation of symmetries of the ABS equations as differential-difference equations lfT2] [TSl fol- 
lows that they correspond to particular cases of Yamilov's discrete Krichever-Novikov (YdKN) equation, or equa- 
tion V4 with V = in Yamilov's terminology ll35l . This equation belongs to the class of Volterra type equations 
and it possesses a local master symmetry |[35l which, in this context, is identified with the first master symmetry 
of the corresponding ABS equation. Here, we present hierarchies of such differential-difference equations and 
their corresponding hierarchies of master symmetries. Moreover, it is shown that the former possess an isospec- 
tral zero curvature representation while the latter possess the corresponding non-isospectral one ||9] [15] and these 
representations are explicitly constructed. 

The organization of the paper is as follows. The next section introduces our notation and reviews some basic 
properties of the ABS equations and their symmetry analysis. Section|3]deals with the construction of conservation 
laws using the first master symmetry while Section 5] contains the higher order master symmetries for equations 
H1-H3 and Q1-Q3. The derivation of the sub-algebra of symmetries in the n direction for equations HI and QIq 
is the subject of Section |5] while the symmetry algebra of HI is studied in Section |6l The relation of symmetries 
to YdKN equation is explored in Section |7] where also isospectral and non-isospectral hierarchies are explicitly 
constructed for these particular cases of YdKN equation. We conclude with an overall evaluation of the results 
obtained in the main body of the paper and several perspectives are discussed. In the Appendix, the ABS equations 
are given along with the basic polynomials related to them. 

2 The Adler-Bobenko-Suris equations 

In this section we introduce our notation and present some properties of the ABS equations. Finally, we review 
some useful results on the symmetries and conservation laws of the latter equations. Further details about the ABS 
equations, their derivation and properties can be found in [2. 3|, their symmetries in ifTT] [TSl l24] l28l |31 J and their 
conservation laws in |[T7ll23l . 

In what follows we will deal with scalar partial difference equations which involve the values of a function 
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ti : H' K, of two independent variables n and m. The dependence of u on n and m will be denoted by 

Un+i,rn+j ■= u{n + i,m + j) , i, j e 

In certain cases and in order to avoid the use of many indices, functions which depend explicitly or implicitly on u 
and its shifts, i.e. functions of the form G{n, m, Un+ii,m+ji, ■ ■ ■ , Un+i^^m+ji), will be denoted simply by a letter 
without indices. In any event, it will be clear from the text if a symbol with no indices depends on the independent 
variables. 

Also, we will use the two shifts operators and the two difference operators corresponding to the two lattice 
directions. The former will be denoted by oS^„ and 5^rn, respectively, while the latter are defined as 

An := - 1 and A™ := - 1 , 

respectively. Finally, the variational derivative of a function g{un+ii^m.+ii, ■ ■ ■ , Un+i^.,m+ji) is defined as 

- - 2^ l^ -^n Q ■ 

p = iiq = jl ^' ^ 

The ABS equations lUO, to which we are interested in, involve the values of u at four different points which 
form an elementary quadrilateral on the lattice and depend explicitly on two lattice parameters a and /3, i.e. they 
are equations of the form 

Qn,m • — Q(^n,m! ''^n+ljjrn ''^n,m,+l ) ^n+l,m+l 1 Q;, /3) — 0. (1) 

The function Qn,m is an affine linear polynomial with respect to its first four arguments and equation ([T]) is invariant 
under the mutual interchange of Un+i.m with Un,m+i accompanied with the interchange of the lattice parameters : 

Qi,Ufi,rm ^n,m+l) ^n+l,m) ''^n+l,m+l) Pi Ck) — i Q{Un,im '^n+l,mj Un^m+lj '^Ti+l/m+lj /^) • 

In terms of the polynomial Qn,m and the discriminant-type operations 

^x,yg{x,y) := gdxdyg - {d^g) {dyg) , &^g{x,y) := {d^gf - 2gd'^g, 
one can define four functions. Namely, 

K{a,/3) f{Un,m,U„+l^m,a) = K,{a,/3)fn,m ■= ^n„,^+i,n„+i,^+i Qn,m , (2a) 
h{'Un+l,rrn 'Un,m+lT 01, /3) •— ^Un.m,Un+i,m+iQn,m i (2b) 
^{'^n.m) = fn,m • — ^Un+i,m fn.m • (2c) 

Here, K{a, /3) is an antisymmetric function of the lattice parameters, fn,m and h are symmetric and biquadratic 
polynomials of the indicated values of u and r„ m is, in general, a quartic polynomial. To make the paper self- 
contained, the ABS equations along with their corresponding polynomials are given in the Appendix. 
All of the ABS equations possess a generalized symmetry [28 J of the form 

^1 ~ ( ~ ~ 9^M,i + l,m/n,m. ) ^Un,m ~- ^«n,m ) (3) 

\Wn,m ^ J 

where 

'UJn,m ■ — ^n+l,m ^n— l,m • (4) 

Additionally, they possess an extended generalized symmetry |[28l 

Gi = nRi du„^^ + ^i{a)da ■ (5) 
The function appearing in Gi is given in the following table. 
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H1,H2 H3 


Ql, Q2 Q3 


Q4 


6(a) 


-1 - 

^ 2 


-1 - 

^ 2 


1^/4:0^ -920-93 



Remark The form of function depends on the particular equation and its parametrization. Here we use the 
parametrization given in ||2] but always one can choose a parametrization in which = 1. In particular for Q4, 
the parametrization used here is due to Adler lH] but one can choose to work in Hietarinta's parametrization iflTI 
at which = 1. □ 

One can construct a hierarchy of symmetries [31] by employing the fact that Gi is a master symmetry of Vi 
Il24l . i.e. [Gi, Vi] / and [[^1,^1],^!] = i]. The resulting hierarchy has the following form lE] 

Vk := [Gi,Vk.i] = Rk ^ (Ri) , k = 2,3,..., (6a) 

where 

00 

^ := ^ p {^PR,) + ei(a) . (6b) 

p=— 00 

From the above definition it is apparent that each Rk depends on {un-k,m, Un-k+i,m, • • • j Un+k-i,m, Un+k,m)- 

Due to the covariance of the ABS equations, a similar hierarchy of symmetries Vi and a master symmetry 
Gi follows from the interchanges {n,Un+i,m+j , C(, (3) — y {m,Un+j^rn+i, (^,0!). It was shown in ||3T1 that the 
members of each hierarchy commute with each other, i.e. 

[Vi,Vj] = [Vi,Vj] = 0, Vi, j G F. 

Moreover, the commutators between members of the two hierarchies and Gi, Gi lead to trivial symmetries, i.e. 
their characteristics are identically zero on solutions of the corresponding equation Qn,m = [22], implying that 

[Vi,Vj] = [Gi,Vj] = [Gi,Vj] = [Gi,Gi] = on solutions of Q„,„ = , Vi, j G M. 

In what follows, we will always assume that commutators of symmetries are evaluated on solutions of Qn,m = 
without writing it out explicitly. 

Another way to construct the above hierarchies of symmetries is to use the pseudo-difference recursion operator 
given in ifTTI . From the latter, one may also construct an infinite sequence of canonical conservation laws by 
considering the residues of this operator. The first conservation law, following from the logarithmic residue of this 
pseudo-difference operator, is 

, /, ^n,m \ .A h{Un-l,m,Un,m+l,a, ^)\ 

Am log „ ] = log / • (V) 

The form of a conservation law implies that one can introduce, at least locally, a potential function and represent 
the conservation law as the compatibility condition of a first order system. For the above conservation law, we 
introduce the potential 4>^^^ = 4>n!ra through the system of equations 

An0(°) = P% := log ^ , A„0(O) = 4% := log ^(""-;--';».-+^'«'/^) , (8) 

/n-l,m K-KCH, p)jn-l,m 

the compatibility condition of which is relation (|7]l. In what follows, we will refer to system ([8j and other similar 
systems, like ([141) and ([TSl l below, as conservation laws as well. 
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3 Conservation laws and master symmetries 

It is well known Q El llll l22l that one can employ symmetries and master symmetries in the construction of 
conservation laws. In this section we follow an approach similar to this in [7 1 which makes use of master symmetry 
Gi and first canonical conservation law ([8]l in terms of which an infinite sequence of conservation laws for all of 
the ABS equations is derived. 

We start with the system constituted by the ABS equation ([T|l and system ^ for the potential (/)(°). This system 
is consistent and its symmetries obviously follow from the corresponding ones of the ABS equation by extending 
them to the (/)('^) -direction. This means that a generalized symmetry of this system has the form 

Vi = Ridu„,^ + <^Pd^,o), (9) 

where ^^^^ is a function to be determined so that 

A„<I>f ) = yniRi)du^^,,^P^^l + Ridu„,^P^^l + ynHRi)du^-^,rnA , (10a) 

Using the identity 

fn,m lo f fn—l,m _, ^ o f 

which holds for any symmetric and bi-quadratic polynomial /, it is a straightforward computation to show that 
the right hand side of equation (llOal ) can be written as A„ (^Ridu„ ^p^-i rnj > which in turn implies that <I>^'^^ = 

Ridun mPn-1 m'^'^rn, whcrc Cm is a function of m. Finally, in view of the latter relation and using equation ([T]), the 
second of equations (fTOl) implies that Cm = cq, where cq is a constant. Since this parameter reflects the translation 
invariance cj)'^'^'^ i— )• (f)^^'^ + e of system ([8]l, we set cq = and, hence, a symmetry generator of system dl}, ([8} has 
the form 

Vl = Rldu„^^ + Rldu„,mPnll,nidc^(0) ■ (12) 



Applying the same procedure with master symmetry Gi by extending it to the direction of potential 

Gi = nRA^^^ + (nRA,,y^l,^^ + ^(i)) d^,o, + 6(a)5, , (13) 
it follows that (p^^^ must obey the system of equations 

^n^^'^ = pffm ■■= ^ (A) - ^n$f ^ , A^,/.(l) = a'^l := ^(4%) . (14) 



In particular, using the identity 

I /• .„ ;) f 

0. (15) 



fn,m du„^fn 



m 



which holds for the polynomials /„ „ corresponding to the ABS equations, function pn^m can be written also as 

(1) _ ^ 

Pn,m ~ '^Un,m^l ■ 

Wn—l.m 
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The compatibility condition of equations ([141 ). i.e. AmPn ln = An<^n}n, is satisfied in view of the corresponding 
discrete equation. Thus, it is another conservation law of the ABS equation Qn,m = 0, which actually coincides 
with the second canonical conservation law given in ifTTll . 

Now we can include (p^^^ in our derivation by considering the system 

Qn,m = 0, A„</.«=/5«^, A^0«=agL, z = 0,l. 

To find its symmetries we have to extend Vi and Gi in the direction of the new potential. These considerations 
result to 

lion (1) I du„„,rn,m '^Un,m'''n,m\ 

+ [Rid..,^Pn-.,r. + - + j d,a^ , (16) 

and 

G, = nVi+ </.(^)a^(o) + + ei(a)9„ , (17) 

respectively, where r„ m is the discriminant of fn,m as a polynomial of Un+i,m, cf. relation (l2c] l. Function (p'^'^'f 
appearing in Gi obeys the system of equations 

where is the -component of Vi. The compatibility condition Ampi^L = A^cri^m of the above system 
is satisfied by taking into account only equation Qn,m = and its shifts, i.e. it is another conservation law of 
Qn,m = 0. Actually, it is equivalen{3 to the third conservation law derived in [flT]. 

We can continue this procedure by including (^(^^ in our system and extending the two symmetry generators 
Vi and Gi in this direction. This leads to the introduction of a new potential ^'^'^^ thi^ough the system 



(2) 

where ^\ is the 9 i (2) -component of Vi, i.e. 



<Pl - iilOu„,^Pn-l,ra + \ 1^2 77., \ + 



Again, the compatibility condition A^Pn^L = A„cji% is another conservation law for equation ([T}- 

We can repeat infinitely many times this procedure and construct an infinite hierarchy of conservation laws 
along with the corresponding functions . In particular, the conserved densities and fluxes are given recursively 
by the relations 

Pn.m =^ (^Pn,m ''^ ~ -^rMx ^ i (^njn ~ ^ {^\i,m ''^ ) ^ = Ij 2, 3, ■ ■ • , (20) 

or, equivalently, 

pSL = (A) - Y.^' (^f "'"'0) ' = (^S^) ' ^ = 1' 2, 3, . . . . (21) 

i=o 



'Two conserved densities are called equivalent if the variational derivative of their difference is identically zero 0171 . 
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(k) 

However, it is not easy to derive a simple recursive formula for functions ^\ . Here we list the first five which we 
were able to find explicitly. 

CF,(0) - VP, JO) 

(1) _ p rj (1) fn,m du„ „^rn^rn , ^Un.m^n 



•Pi - ^lC«n,™/'n-l,m + ^-2 + 24 

;(2) _ Tj fl (2) \ (S> ( ''^n,m du„„^rn^m , '^Un.m^r, 



$1 - i^ld„„,„,p„_i_^ + =^ —2 — — + 



2w^n-l,m AWn-l,m 24 



1>i - iiid„„_,„p„_i_^ + ^ —2 — — — + 





,m 




1'n,m 


2^1 


-l,m 


4u;„- 


-l,m 




,m 






2^1 


~l,m 


AWn- 


-l,m 



24 y ^n-l,m^n-2,m 

Only for two equations, namely HI and QIq, for which r„ ^ = 0, one can generate recursively the whole hierarchy 
of ^f^^ and this construction will be given in Section|5] 

(k) / N (k) 

It is obvious that the conserved density pn,m depends only on {un-k-i,mT ■ ■ ,Un^m,Un+i,m.), while ^\ 
depends on {un-k-2,m 1 Un—k—i,im • • • 1 Un^m^ ^n+i,m)- We Can provc inductivcly that this construction does not 
terminate. The first five conserved densities are neither trivial nor equivalent to each other, which can be checked 

(k—l) 

directly using the above formulas and any of the ABS equations. Suppose that after k iterations, density pn^m is 

(k) (E) 

neither trivial nor equivalent to any other lower order density, but the next one pn,rn is equivalent to pn'm, where 
i < k. This means that 6u i^Pnjn — Pn,m^ = 0. But the only term in the difference of these two densities which 

involves the value Un-k~i,m is —k.y^^ (Ri) du„_^ ^Pn,rn'^ 7^ and its contribution to the variational deriva- 
tive is non-zero, since „ (Rl) 9„„„,,„pSm^^)) = /n+l,m (^"n+i.^/'iVfc+l.m) /^n+l,m / 0- 
(k) (k) 

Thus, pn,m is not equivalent to any lower order density. From the same arguments, it follows that pn,m cannot be 
trivial as well. 

4 Hierarchies of master symmetries 

Once we have constructed the hierarchy of potentials (f)^''^ and the symmetry generators 

k 

Gi = ni?i5„„,„ + ) + c/^^^+D) d^^,, + , (23) 

k 

we can find how the higher order symmetries Vn ^ act on the potentials (/i^^). This can be done by considering the 
commutators of Gi with these symmetries starting with Vi. To this end, the symmetry hierarchy is given by 

Vt = Rl dur.,m + ' ^ = 1, 2, 3, . . . , (24a) 

k 

where 

<I>f ) = ( ^ T M {-lY+^-^M^ (^<i>f +^-^-1)) , £ = 1, 2, 3, . . . , A: = 0, 1, 2, . . . . (24b) 
i=o V / 
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Using the hierarchy of potentials cj)^ ' , one can derive a corresponding hierarchy of non-local master symme- 
tries for equations H1-H3 and Q1-Q3. More precisely, this can be done starting with Gi above and 

G2 = (nR2 - 4>^'^Ri - -^Jb^\ + 6(a) 5a + C2{a,(3)dp , (25) 

where functions 6(a), C2 (a,/?) appearing in G2 are given in the following list. 

HI : (a-/3)(C2-6) - 1 = =^ 6 = 0, 6 = ^^ 

a — p 

H2 : 6 = 0, C ^ 



a- 13 



H3o : 4(a2 - f3^){aC2 - /36) - al3{a^ + P^) = =^ 6 = "t- C2 



4' 2(a2 -/32) 



■tl-aiS • 42 — 7 , C2 



4 ' 2(a2 _ /32) 



Qlo : aC2-/3(2-p(- + ^—)=0 =^ 6 = --, (2 = ^^-- 

\a a — pj a a — p a 

QI5 ■ ^2 = , C2 = 

a a — p a 

Q2 : 6 = --, C2 ^ ^ 



Q3 : 6 = - , , , , C: 



a a — (3 a 

a{a^ + 1) ^ a2/3(/32 - 1) 



4(a2_i)' 2(a2 - l)(a2 -/32) 

Remark In the case of HI, H3o and QIq, functions 6> 6 ^re not completely determined and one has to choose a 
solution of the corresponding linear algebraic equations. In the list, a particular solution for each of these equations 
is given and this solution will be used in what follows. Our choices for H3o and QIq coincide with the correspond- 
ing functions of their (^-counterparts, while for HI with the ones for H2. □ 



Remark on the derivation of G2 For HI the characteristic of G2 was derived from the first master symmetry Gi 
by applying to the latter the recursion operator given in [17 |. Then the linear algebraic equation for functions 6 
and 6 follows from the requirement the resulting function to be a symmetry of HI. The same approach was used 
for the remaining equations admitting first order recursion operators, namely equations H2, H3, QI and Q3o. For 
Q2 and Q35, the form of G2 was verified by direct calculations. □ 

One can verify that G2 is another symmetry of the corresponding ABS equation in a straightforward manner. 
Specifically, it is sufficient to show that the infinitesimal criterion 

E E f^^2 - ^^'^Rl - -^ii^) ^^"'"^ + 65.Qn,™ + C2dpQn,^ = 

holds on solutions of Qn,m = and its shifts, and by taking into account that 

^ ^(1) _ Ji) ^ ^(1) _ ^(1) 
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Moreover, the commutator of G2 with any of the symmetries Vj will always produce local symmetries. Indeed, 
writing out explicitly the commutator [G2, Vj] and using the facts that all symmetries commute [Fj, V^] = and 
An'/'''^^ = Pn}m, One finds that the characteristic of the resulting generator has the following form. 



k=-] ^ ' k=-] 

^ (^) -E (e^^Sl) ^+ E (e^>S 

k=-2 ^ n L,m/ y.^Q y n+«,m y.^^ 



du. 



n+k,m 



All the functions involved in the above expression are local, i.e. do not depend on the potentials 4>^'^\ implying that 
indeed [G2, Vj] is local. Moreover, it does not depend explicitly on n. 

Additionally, G2 is another master symmetry for Vi since it can be shown that 

for HI, H2, Ql and Q2 : [G2, ^i] = ^^3 , 

for H3 and Q3 : [G2, Vi] = + ^Vi , 

implying in each case that [[G2, Vi\,Vi\ = 0. Finally, it can be verified that the action of G2 on conservation law 
O leads to the second conservation law (fTSl ). i.e. G2 produces higher canonical conservation laws in the same 
fashion as Gi does. 

One can construct the next master symmetry by considering the commutator of Gi and G2. At this point it is 
more convenient to consider separately equations HI, H2, Ql and Q2 from equations H3 and Q3. For the former 
equations one finds that the higher order master symmetries have the following forms. 

+ Ma)da,ii-i{a)da, + Q_i{a,l3)dp], £ = 3,4,5,---. (26) 

For equations H3 and Q3, the resulting symmetries partly follow this pattern. More precisely, the next master 
symmetry G3, following from the commutator [Gi, G2], reads as follows 

Gs = (l:R3 - ^Ri - - <^(')i?i - ^ ) ) du,, ^ + [ ei 9„ , 6 5„ + C2 5^ ] , (27a) 

V2 12 \2Wn-l,mJ J 



while the next two symmetries G4 and G5 follow the general pattern (I26l l. Continuing with the next two symme- 
tries, one finds that 

( - - ( ^ ) C^^'-'^R, + l - { 2Wn-lm ) J ^ i^^^^^ ^^^a + C^dp] (27b) 



and 
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respectively. 

It can be shown by induction that the commutators [Gj, Vj] result to local symmetries. As we have shown, this 
is the case when i = 1 and i = 2. Using the definition of G3 and the Jacobi identity, the commutator [G3, Vj] can 
be written as [Gi, [G2, Vj]] — [G2, Vj+i]- Since Gi is local and we have shown that [G2, Vj] is local for any i, it 
follows that [G3, Vj] is local as well. Assuming that [Gk, Vj] is local for some A; > 3 and for all j and using the 
same arguments, it follows that [G^+i, V,] = [Gi, [Gk, Vj]] — [Gk,Vj+i] is also local. 

Remark In the continuous case there exist sufficient conditions which guarantee the locality of symmetries re- 
sulting from the commutator of master symmetries with other symmetries [26|. More precisely. Proposition 2 in 
||26]| gives the necessary conditions on the form of the non-local part of a master symmetry so that to produce local 
symmetries. Actually, comparing the non-local part of all the above master symmetries with the form assumed in 
Proposition 2 of [26] . it is obvious that G^'s are of this form, with Wy and uj^ in [26] identified respectively with 
and (/)(^) here. □ 

Once the symmetry sub-algebra n generated by Vi and G^ has been constructed, another symmetry sub-algebra 
m follows from the covariance of the ABS equations. Specifically, hierarchies of symmetries V, master symmetries 
Gi and potentials follow from the usual interchanges {n,Un+i,m+j-,ot, (3) — > {m,Un+j,m+i, P,Oi) which 
generate the symmetry sub-algebra m. These two sub-algebras, along with the point symmetries each equation 
possesses, form the symmetry algebra t) of the corresponding equation. To find the commutators among the 
elements of the two sub-algebras, one has to extend the action of Vk and Gk in the direction of the potentials c^*^*) 
and then evaluate the commutators. However, to find general expressions valid for any equation in this class is 
very difficult. This can be done for each equation separately and for equation H 1 the corresponding results are 
presented in Section |6] 



5 Discrete potential and Schwarzian KdV equations 

Having constructed a sequence of symmetries and master symmetries in the manner described in the previous 
sections, one interesting question is what is the algebraic structure of these symmetries. In the case of differential 
and differential-difference equations usually one finds that all these form centerless Virasoro type algebras. In this 
section we show that the symmetries (V^, Gk) of HI, as well as of Ql with 6 = 0, form a centerless Virasoro 
type algebra. This is achieved by finding explicit expressions for all symmetry generators and their action on the 
potentials (f)^''\ and, consequently, computing their commutators. Due to the covariance of HI and QIq, similar 
relations hold for the symmetries (Vk, Gk) in the m direction. 
The lattice potential KdV equation (also known as HI) 

{Un,m - Un+l,m+l){Un+l,m - Un,m+l) - « + /3 = (28) 

and the lattice Schwarzian KdV equation (also referred to as QIq) 

a(^n,m — Un,m+l){Un+l,m. — Un+l,m+l) — l3{Un,m — Un+l,m){Un,m+l — Un+l,m+l) = (29) 

are the only equations in the ABS classification lO for which the polynomial Vn^m is identically zero. For these 
two equations, it is easier to construct the hierarchies of conservation laws and hierarchies of symmetries using the 
results of the previous sections. 

The characteristic Ri of the first symmetry Vi is 

for m: Ri = ^— , (30a) 
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and 

lOr l^lg : Ki = , (jUDj 

respectively. Now, starting with the conservation law (H) which takes the following form 

£ TT1 fO) 01 \ fO) 1 ('Un— l,m ~ Uji rn+l) ^ 

for HI : p'^"l^ = 2 log [un+i.m - Un-i,m) , cTiiL = log '—^ ■ , (31a) 

p — a 

and 

forQlo : p'^^>^ = 2\og ^ + log a, cj),^/„ = 2 log + log^^— ' (31b) 

(k) 

respectively, one can construct the hierarchy of conserved densities and fluxes along with functions . Specifi- 
cally, relations (l20l ). which define recursively densities and fluxes, now simplify to 

(fe) ^ / (0) \ (k) ^ ^k / (0) \ .32a) 

where is the "negative" part of i.e. 

-1 

J2 ^^i(^i)9„„+,,^ - (32b) 

£=-oo 

In terms of the above densities, functions <I>^ can be expressed as 

= Ridu^,^P^^i,^ , i = 0, 1, 2, . . . . (32c) 
The hierarchy of symmetries is given by 

Vk = Rk du^_^ + ' A; = 1, 2, 3, ... , (33a) 

i 

where are defined in (13 2c I ) and, in terms of which, <I>^*^ are determined by relation (124b l i which now becomes 

^k^ = E(^~,') i-^r^-'^' ''0 • ^^^'^ 

i=o V / 

On the other hand, the hierarchy of master symmetries is given by 



and 



+ 6(a)c'„ + Cfc(a,/3)5^, k = 2,3,..., (34b) 
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with 



{a-l5Y 



for HI : efc = 0, Cfc = ..fcl^ , (34c) 



(35) 



-(fc-2)! , (A; -2)! (A: - 1)! 

a'^ ^ (a — pj"^ ^ a'^ ^ 

Computing the commutators of the above symmetries, one is led to the following relations. 

[Gi^Vi] = yj+i , [Gk, Vj] = ^^Vk+j , A: = 2, 3, ... , 
[GuGj] = , [G„ Gj] = (j - , i,j >1. 

These suggest to rescale the symmetry generators as 

% = jkkjlVk, k = l,2,..., 
'^1 = Gi, % = ik^Gk , k>l, 

in terms of which the above commutation relations become 

= j^+i, (36) 

These relations imply that the algebra n of the above symmetry generators is the semi-direct sum of the abelian 
algebra {^ j and the centerless Virasoro algebra {'^k}, n = {^} (Bs {%}■ As consequence of the covariance of 
equations HI and QIq, the symmetry generators will satisfy the same commutation relations and, hence, 

will form the same algebra. 

Remark Working out explicitly some of the commutators for equations H2, Ql^ and Q2, it is apparent that their 
symmetries satisfy similar commutation relations suggesting that they form, after a suitable rescaling, similar al- 
gebras. However general expressions, like (I32ti34l i. are difficult to be derived in terms of which this statement can 
be proven. □ 



6 The symmetry algebra of HI 

The symmetry algebra f) of HI is generated by its point and generalized symmetries. The first ones have been 
found in [24, 28] and the generalized symmetries were the subject of the previous section. Here, we will explicitly 
compute the commutators among the generators of these three sub-algebras and, consequently, find the structure 
constants of f). 

Let us begin with the sub-algebra p of point symmetries generated by the vector fields ll24] |281 

j>i i>i 
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One easily verifies that the only non- vanishing commutators among them are the following ones. 



[^1, >0] = =^2 , = , [%, ^l] = \^1. ^2] = \^2- 

Remark The reason for denoting two of the above generators by and % is that these two specific symmetries 
are related to the hierarchies of symmetries "fl, % and master symmetries Wi, Wi, respectively. As it was noted 
in Section m ^2 can be derived from ?fi by applying to the latter the recursion operator which was given in IfTTl . 
namely operator + l)A^^i?i(o$^„ — ^n^)- In exactly the same fashion, can be constructed from % by 

applying to it the same recursion operator and then using the arguments mentioned in SectionlH Hence, the scaling 
generator % is the seed symmetry from which the whole hierarchy of master symmetries is derived by applying 
successively to it the pseudo-difference recursion operator + l)A~^i?i(o5^„ — ^^^). Similarly, the result 

of the application of the same operator on Yq is the first symmetry Yi implying that the former can be taken as the 
seed symmetry of the hierarchy □ 

Before considering the commutators among generators of the three sub-algebras, let us rescale the generalized 
symmetries and re-define them as follows. 

Wk := (a - /?)^ % , M := (a - % , ..7. 

It should be noted that since a and j3 are parameters, we can freely multiply any symmetry by a function of 
the lattice parameters with the resulting vector field being again symmetry generator of HI. For instance, the 
infinitesimal criterion for Mk to be a symmetry of HI is the same with the one for differing by an overall 
multiplicative factor {a—j3)^. The commutation relations among the generators Wk, Mk of n and the corresponding 
ones of m follow from (l36l ). (|37] ) and read as follows. 

= j - ^i) = j {^i+j - ^j) (38) 

It is straightforward to verify that all the point symmetries commut^ with the above generalized symmetries, 
i.e. [p,n] = [p,m] = 0. Thus, symmetry algebra \) of HI is generated by the sub-algebra p of point symmetries 
and the sub-algebra g of generalized symmetries, and, in particular, = p © g. For the algebra g, we know 
that g = n + m as a vector space. Now, we will try to find the commutators among the generators of these two 
sub-algebras. 

For this purpose, we will extend the action of "Vx, 1^, and in the direction of the potential ^^^^ . Once 

these actions are found, the covariance of HI implies the corresponding actions of >i, and on i;^^^). 

Finally, these results combined with definitions (1371 ) will deliver the corresponding action of Wi, "Wi, M'i and M'i 
for i = 1,2. 

The action of 1^ and on (^^^^ is given by 

; (1) ^ _ 1 / 2 / 1 '"n.m-l - ^n-l,m ^ 1 ^ o 



^Although all the other point symmetries commute with the hierarchies of '/i and 'i^i, scaling symmetry Sfo satisfies [l^o, "^j] = j ^ 



and [Wo/'^j] = j • From the latter and definitions J37) follows that commutes with Wi and .Jifi. 
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and n i^^'' , respectively. In the above relation, Wn,m. is the function defined in ^ and 

Vn,m ■ — 'Un,m+1 'Un,m—1 ■ 

The action of ^ on (f>^^^ can be derived from the definition of the latter, ^2 = ['^i, '^i], and it is 



Finally, ^2 acts on (^(^^ as 



- ' - ) 

Now, the commutation relations 

[J^ , = J^i — J^i , ['^^ , = 2 — t;^^ (39a) 

and 

[^2, #1] = -n, [^2, ^2] = 4 (j^2 - ^2) + 3 (jTi - ^1) (39b) 
can easily be verified. From commutation relations (l38T l. (|39l) and the Jacobi identity it follows that 

= 0, 

L^'d = + E ^M^jj ' (40a) 

i=i ' i=i 



a. - ( ) • (40b) 



where 

i + i 

The above commutators hold for any fc, ^ = 1, 2, 3, . . ., and along with the ones in (l38l) describe the sub-algebra g 
of the generalized symmetries of equation HI. 



7 Symmetries and differential - difference equations 

The relation among symmetries of the ABS equations and differential-difference equations of Volterra type is well 
known |[T2l [181. In particular, the former can be regarded as particular cases of Yamilov's discrete Krichever- 
Novikov (YdKN) equation f34l[35l . It is also known that YdKN equation admits a master symmetry in terms of 
which the whole hierarchy of symmetries can be constructed ll27l . In this section, we present for the particular cases 
of the YdKN equation corresponding to symmetries of H1-H3 and Q1-Q3, the hierarchy of master symmetries 
along with their Lax pairs and discuss the algebraic properties of latter [i9iil5 J. 
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In this interpretation, we assume that u also depends on an infinite number of time variables tk and r^, which 
essentially correspond to the group parameters of the transformations generated by the symmetry generators 
and Gk, respectively. The symmetry condition now translates to the requirement the flows 

^ = Rk (41) 

dtk 

to commute with the discrete equations under consideration, i.e. dt^Qn,m = 0. The latter relation holds on 
solutions of Qn,m = by taking into account (|4T]) . Moreover, from the construction in Section [3j it follows that 
the relations 

dt.pSm = ^n^t^ , A: = 1,2,..., i = 0,l,..., (42) 

hold on solutions of (|4T1) . In other words, they are conservation laws of the corresponding differential-difference 
equations (|4TI ). Actually, for A; = 1 and up to rescaling, the first two conservation laws, which correspond to 
i = and i = 1, respectively, coincide with the first two integrability conditions appearing in the classification of 
integrable Volterra type equations, cf. relations (A.ll) in [4|. On the other hand, the third conservation law, i.e. 
for i = 2, is equivalent to the third integrability condition given in [4]. 

In the same fashion, for master symmetries Gk corresponding to H1-H2 and Q1-Q2, we have that the flows 



dUn,: 



dri 



n Ri , (43a) 

-2 



accompanied with 



(43b) 



|^ = efc(«), |^=a(«,/3), k = l,2,..., (43c) 
dTk dTk 

commute with the corresponding discrete equation, i.e. dr^Qn,m = holds on solutions of Qn,m = in view 
of (l43l) . For equations H3 and Q3 similar considerations hold using the corresponding characteristics of master 
symmetries (|27] ). 

Concentrating on the hierarchies of differential-difference equations (|4T]) and their master symmetries (l43T l. 
we will see how the former possess an isospectral zero curvature representation while the latter possess the corre- 
sponding non-isospectral one fO'.T?!. In other words, we will construct Lax pairs for all these equations in which 
the spectral parameter, the role of which is played now by the lattice parameter /3, will be considered as a function 
of the time variables Tk and independent of tk- Specifically, the second relation of (I43cl ) will provide us with the 
particular dependence of the spectral parameter /3 on the time variables r^. 

To start with, let us construct the spectral problem to which these hierarchies are related to. It has the form 

'^n+l,m = L ^n,m , (44a) 

where 

_ \ f t^ttn + l,m + lQn,ni ~du„ + l,m + l^U„^rn + lQn,m J (44b) 



K{a, I3)fn,m \ Qn,m du„ ,^j^^Qn;, 

Here, Qn,m and its derivatives are evaluated at n„ m+i = Un+i,m+i = and j3 plays the role of the spectral 
parameter Matrix L is coming from the corresponding Lax pair of the ABS equation to which equations (|4T]) . (I43l) 
are related to. In particular, matrix L follows from the multidimensional consistency of the corresponding ABS 
equation Qn,m = Il6l|20l and the above form of matrix L was given in OTi 
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The isospectral hierarchy can be derived systematically and has the following form OT 



dtk 

The matrices M^. are given recursively by 



Mfe^'.^, A; = 1,2,.... (45a) 



with 

Ml 



Mk=M''~^Mi, A; = 1,2,..., (45b) 
( 77^^ 7;du„+i,m] (-/n,mL-^a„„_„L) . (45c) 

\Wn,m ^ J 

The compatibility condition + LM^ = (^„Mfc) L of (l44b . (1451 ) is satisfied in view of equation (|4T]) . 

On the other hand, the non-isospectral hierarchy follows the pattern of master symmetries with being 
replaced by matrix and the last term involving polynomial r„ m being substituted by a suitable matrix H. More 
precisely, the non-isospectral hierarchy corresponding to (|43] ) is given by 

= Ni^„,„ := (nMi) (46a) 

along with 



Wn^m := f - 5^ ( ^ ^. ^ ) 0(^-^-i)M,+i -^^-2 (H) 1 ^ > 2, (46b) 



|^=eKa), |^ = CK«,/?), ^ = 1,2,3,.... (46c) 



On the other hand, for H3 and Q3 in particular, we have 



5t6 



9r7 





— M 
12 


(^-M„ - 


n ^ ^ 
— M. 
30 ' 


-Mr - 

U 


, ^ 
— M- 
12 ' 



Mi-^Hj*„,„, (47) 
-j)M,+i-^4(H)j (48) 

-^■)M,+i-i^5(H) I (49) 



Matrix H appearing in the above relations can be derived, up to a point, in a systematic way and its final form 
depends on the particular equation. These precise forms are given in the following list. 

HI, H2, Ql : H := p^^l^Mi - 2 (^'^iJi) a„„_,,„Mi - d^Mi (50a) 

H3, Q3 : H := p^^lu, - 2 (^-^iJi) a„„_,,„Mi + | a„Mi + ( ° ? ) (50b) 

Q2 : H := p^^M, - 2 {y-'Ri) 9„_,„„Mi - 9„Mi + ° ) (50c) 

The compatibility condition dr^.L + LN^ = {^n^k) L of (l44l ). (|46] ) is satisfied in view of equation (l43l) . 
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Remark We would like to emphasize that the formulae for the hierarchies of differential-difference equations 
(|4TI ). their canonical conservation laws (02), as well as their Lax pairs (|44] |. (l45T i. are valid for the symmetries 
corresponding to any of the ABS equations, i.e. including Q4 as well. However, only the formulae for the first 
master symmetry (I43al i and its Lax pair (I46al i are also valid for Q4 [i311 . □ 



As explained above, matrices L, Mi , Ni and N2 were derived in a systematic way. The remaining matrices 
follow from the consistency of the different linear equations for the potential ^'„,m- More precisely, the compat- 
ibility condition of equation (l45l) with k = I, and equation (I46al ) yields M2, while the compatibility condition of 
(|45] ) with k = 2, and (I46al) is matrix M3. Thus, the matrices M^ of the isospectral hierarchy are constructed in the 
same way hierarchy (1411 was derived. 

Mfc+i = [Mfc,Nil := 1^ - ^ + [Mfc,Ni] , k = l,2,... . 

Similarly, the non-isospectral hierarchy can be generated in the same manner starting with equations (I46al l and 
(I46bl ) with i = 2 and their consistency. This approach leads to 

Nfc+i = [Nfc,NiI := ^-^ + [Nfc,Ni] , A; = 2,3,... . 
on dTk 

Finally, the commutativity of hierarchy (1411 implies that 



IM,,M,1 :=^-^ + [M.,M,]=0. 

Oti Oth 



In all the above relations the derivative with respect to r^j is understood as 

d X; ^^^1 +Ua)d^ + Q,{a,P)dp, k = l,2,.... 



The bracket |-, •]] was introduced in lITSl where the general algebraic structure related to discrete zero curvature 
equations was developed and was illustrated by the examples of Volterra and Toda lattice hierarchies. In particular, 
in those examples it was shown that the Lax operators form Virasoro type algebras. In fact, employing the analysis 
and the results of Section |5j it can be verified that the following relations hold for the hierarchies of Lax operators 
corresponding to equations HI and QIq. 

[M„MJ = 0, lMj,Nil = jM,+,-, [N,-,Nil = (j-i)Ni+,. (51) 

For equations H2, H3 and Ql^-QB, we remark that one could employ the results of lITSl and use the fact Il35l 
that, under some easily verified conditions, the first member of hierarchy (1411 is related through two different 
Miura transformations to either the Volterra equation dtUn = Un{un+i — Un-i), cf. also [12], or the Toda system 

dtUn = Un{Vn+l - Vn), dtVn = Un - Un-1- 



8 Conclusions and perspectives 

We have presented a constructive approach to the derivation of symmetries and canonical conservation laws for 
the integrable difference equations of the ABS classification. The basic ingredients in this construction are the 

first symmetry Vi, the first master symmetry Gi and the first canonical conservation law ipn}n, 0"!% ) in terms of 
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which the rest members in these hierarchies can be derived systematically. The approach for the construction of 
canonical conservation laws can be applied to any equation and it would be interesting to apply this method to the 
non-autonomous deformed H equations [3] the symmetries of which were studied in [33], as well as to the most 
general equation in the class of Kleinian invariant equations, i.e. equation Q4 with arbitrary parameters ||29l . 

For equations H1-H3, Q1-Q3 the second master symmetry G2 was given explicitly and it was used with 
Gi in the construction of a corresponding chain of master symmetries. All the higher order master symmetries, 
including G2, are non local since their characteristics depend explicitly on the potentials </)*^*) corresponding to 
canonical conservation laws. It was proven that these non-local master symmetries map local symmetries to local 
symmetries. In particular, for equations HI and QIq, it was shown that their symmetries and master symmetries in 
each lattice direction form centerless Virasoro type algebras, and the whole symmetry algebra for HI was derived. 
The construction of the second master symmetry of Q4 is still an open question. An interesting question related to 
master symmetries and conservation laws is if these master symmetries can be regarded as variational symmetries 
of the ABS equations, a Lagrangian of which is given in 1,32 J . and if these canonical conservation laws can be 
derived by implementation of Noether's theorem. 

The interpretation of symmetries as differential-difference equations led us to the derivation of hierarchies of 
integrable differential-difference equations along with corresponding hierarchies of master symmetries and canon- 
ical conservation laws. We also presented isospectral and non-isospectral hierarchies to which these differential- 
difference equations and their master symmetries correspond to. It would be interesting to derive similar structures 
for the general YdKN equation, i.e. the first generalized symmetry of Adler's equation, which is not related to 
Volterra or Toda systems fil2ii35j . 
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A The ABS equations and their polynomials 

With the identifications 

the Adler-Bobenko-Suris equations are given in the following list. 

HI {u- z){x -y) - a + ^ = 

H2 {u - z){x - y) + {13 - a){u + X + y + z) - + l]"^ = 
H3 a(ux + yz) - I3{uy + xz) + 6{a^ - ^'^) = 
Ql a{u -y){x - z) - l3{u - x){y - z) + d'^a/3{a - j3) = 
Q2 a{u — y){x — z) — j3{u — x){y — z) + 

a/3{a - P){u + x + y + z) - aP{a - /3)(q^ - a/3 + = 
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Q3 (/3^ - a^)(nz + xy) + /3(a^ - l)(na; + yz) 

-a(/32 - l){uy + xz) - d^ia"^ - /3^){a^ - l){l3^ - 1)/(4q/3) = 

Q4 aQUxyz + ai{uxy + xyz + yzu + zux) + 02 (nz + xy) 

+a2(ttx + yz) + d2{uy + xz) + 03(1* + x + y + z) + a4 = 

In Adler's parametrization lH], the parameters involved in Q4 are determined by the relations 

ao = a + b , ai = —aP — ba , a2 = a0^ + ha^ , 

ab{a + b) 2 (n 2 52 ^, ~ ab{a + b) 2 ^ofl2 92. 

2 4 lo 

with 

= p[a) , 6^ = , p(x) = - 52a; - 53 • 
In Hietarinta's parametrization fTT\, the parameters in Q4 have the following forms 

ao = = k sn{a) sn(/3) sn(a — /3) , oi = 03 = 0, 
02 = sn(/3 — a), 0,2 = sn(a), 02 = — sn(/3), 

where sn denotes the Jacobi elliptic function sn{a\k) with modulus k. 

The polynomials related to the ABS equations are given in the following list. 



Equation HI 
Equation H2 



Equation H3 



• Equation QI 



Equation Q2 



fn,m = l, h{x,y,a,f3) = {x - yf, K{a, f3) = f3 - a, r„,. 



fn,m = '^iUn,m + Un+l,m + «), h{x, y, a, f3) = {x - yf - {a - I3f 
K{a,(3) = (3 - a, rn,m = 4 



fn,m = Un,mUn+i,m + OiS, h{x, y, o, /3) = {ay - /3x){f3y - ax), 

K^a, /3) — a /3 , fn,m, — ^n,m 



fn,m = {{Un,m - Un+l,m)'^ - C?b^)la, 

h{x, y, a, /3) = a/3 {{x - yf -{a- (if 5^) , 
K(a,/3) = -a/3(a - /3), r„,m = 45^ 



fn,m — ((^n,m, "Un+l,™) 2a {Uri,m ~\~ '^n+l,rn) ~\~ Ci )/a, 

h{x, y, a, (3) = a/3 ((x - y)^ - 2(a - ^)\x + y) + {a - /3)^) , 
K{a,/3) = -al3{a - l3), rn,m = ^^Un 
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• Equation Q3 





h{x,y,a,l5) = ^ (4a/3(ay 




• Equation Q4 (Adier's parametrization) 




Polynomial h is omitted because of its length but it can be derived directly from its definition (l2bl ). 
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